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In this work, we investigate the effect of disorder on the topological properties of multichannel 
superconductor nanowires. While the standard expectation is that the spectral gap is closed and 
opened at transitions that change the topological index of the wire, we show that the closing and 
opening of a transport gap can also cause topological transitions, even in the presence of nonzero 
density of states across the transition. Such transport gaps induced by disorder can change the 
topological index, driving a topologically trivial wire into a nontrivial state or vice versa. We 
focus on the Rashba spin-orbit coupled semiconductor nanowires in proximity to a conventional 
superconductor, which is an experimentally relevant system, and obtain analytical formulas for 
topological transitions in these wires, valid for generic realizations of disorder. Full tight-binding 
simulations show excellent agreement with our analytical results without any fitting parameters. 

PACS numbers: 74.78.Na, 74.45.-|-c, 71.23.-k 


I. INTRODUCTION 

Topologically nontrivial phases are exotic states of 
matter that have an electronic band gap in their bulk 
and protected gapless excitations at their boundaries.^i^ 
Superconductors, being quasiparticle insulators, also fea¬ 
ture topological phases with a quasiparticle gap in the 
bulk and excitations at their edges. For ID systems, 
these edge states are fermionic zero-energy modes called 
Majorana states.l^lf^ These states attracted intense at¬ 
tention owing to their non-Abelian nature, which led to 
proposals t o use them as topological qubits immune to 
decoherence.l21i2l Although predicted to appear in exotic 
condensed matter s^tems with unconventional supercon¬ 
ducting pairing,recent proposal^^^^ll^oi involving hy¬ 
brid structures of more conventional materials have ap¬ 
peared.!^ This led to the recent conductance measure¬ 
ments done on a proximity coupled InSb nanowire 
which showed possible evidence of Majorana end states 
in the form of zero bias conductance peaks. Other 
experiments reported further observations of zero bias 
peaks in similar settings.!2^!!^ Very recently, scanning¬ 
tunneling spectroscopy experiments carried out on mag¬ 
netic adatom chains on a conventional superconductor 
reported ZBPs at the ends of the chains.^S While it is 
compelling to interpret the observation of these ZBPs as 
signatures of Majorana states, the issue is still under in¬ 
tense discussion.SS 

Semiconductor nanowire structures that are proximity- 
coupled to superconductors are technologically attractive 
platforms for Majorana physics. However, disorder has 
been prominently present in all such experimental sam¬ 
ples to date. This led to a renewed interest in disordered 
superconducting wires, particularly focusing on the ef¬ 
fects of disorder on Majorana states.l^SlESl These works 
focused mostly on disordered p-wave superconducting 
wires (PW wires) and showed that disorder is detrimen¬ 
tal to the spectral gap as well as to the formation of 
Majorana fermions in both strictly ID system^SSHSl and 


in multichannel wires.^^^lHIIII] jn a recent study on the 
experimentally relevant hybrid structures with Rashba 
spin-orbit interaction (SOI) proximity coupled to an s- 
wave superconductor (RSW nanowires for short), some 
of us showed that disorder need not be detrimental to 
and in fact can even create topological order in strictly 
ID wires.!^ We are not aware of a systematic study of the 
effects of disorder on the phase diagram of multichannel 
RSW nanowires. 

In Majorana experiments, the subband spacing is typ¬ 
ically considerably larger than the Zeeman splitting. For 
example, in InSb nanowir es a s ubband spacing of order 
I5meV has been measurecP^Ml together with a g-factor 
of 40 to 58. Zero bias peaks that might signal Majorana 
fermions in these works ar e typi cally measured at mag¬ 
netic fields from O.lmT - 11*221111 and exceptionally up to 
2.5T. In all of these cases the Zeeman splitting remains 
smaller than the level spacing. Hence, one can argue 
that RSW nanowires are more experimentally relevant 
than PW nanowires, which require Zeeman splitting be 
much larger than level spacing. 

In this Manuscript, we investigate topological prop¬ 
erties of disordered multichannel RSW and PW super¬ 
conductor nanowires. The usual expectation for these 
nanowires is that if their topological state is switched by 
modifying certain external parameters (such as gate po¬ 
tential or magnetic field), the spectral gap will close and 
open concomitantly with this transition. We show that 
for disordered nanowires, the closing and opening of a 
transport gap can cause further topological transitions, 
even in the presence of finite density of states (DOS), 
extending our earlier work on single channel wire^^ to 
multichannel wires. We derive analytical expressions for 
the boundaries of the topological phases of a disordered 
multichanneled RSW nanowire and find new topological 
regions in the phase diagram that show up as additional 
reentrant behavior in the experimentally relevant param¬ 
eter regimes. In particular, new topological regions that 
show up in the low magnetic field limit, requires full de- 
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FIG. 1. The multichanneled nanowire of width W, which is 
an RSW topological superconductor with a Gaussian disorder 
having an average value {V) — 0. a) In the leads, we take 
aso, A and V{x,y) to be zero, making the leads metallic. 
Our analytical results assume a semi-infinite wire {L —>■ oo), 
whereas in our numerical full tight-binding calculations we 
use wires of length L 3> Iso- b) The form of the wire 

used to construct the Majorana solutions in section [II A[ The 
part of the wire left of the scattering region is again metallic. 


scription of all spin bands as shown by our analytical 
results (see Fig. [^. Hence, our results go beyond a sim¬ 
ple p-wave description that requires a fully spin polarized 
wire. Finally we perform numerical simulations using a 
tight-binding (TB) approach and find very good agree¬ 
ment with our analytical formulae. 

This Manuscript is organized as follows: We begin 
the next section by specifying the system in question. 
We then derive the topological index in terms of the 
Lyapunov exponents and the effective superconducting 
length of the disordered multichannel RSW wire in sub¬ 
section |II A[ In subsection |II B[ we analytically calcu¬ 
late this topological index using experimentally relevant 
system and transport parameters and compare our re¬ 
sults with numerical tight-binding simulations. We then 
present our conclusions, finding that in disordered mul¬ 
tichannel RSW nanowires with experimentally relevant 
parameters, the topological phase diagram is fragmented 
and previously unreported reentrant topologically non¬ 
trivial regions appear. In the Appendices, we detail 
the calculation of the mean free path of the system 
(Appendix detail our numerical simulations (Ap¬ 
pendix]^, present a full bandwith versions of our plots in 
the main text as opposed to the low energy region (Ap¬ 
pendix 0. and finally present our plots similar to the 
RSW system but preoduced for a p-wave nanowire with 
disorder, as system previously studied in literature, for 
completeness and comparison (Appendix [Pj). 


II. TOPOLOGICAL ORDER IN DISORDERED 
MULTICHANNEL WIRES 


In this section, we investigate the topological prop¬ 
erties of multichanneled topological superconductor 
nanowires. Such wires are experimentally realized 
by proximity coupling a semiconductor nanowire with 
Rashba spin-orbit interaction to an s-wave superconduc¬ 
tor (RSW, see Fig. (a)). The quasiparticles in RSW 
nanowires are described by the following Bogoliubov-de 


Gennes (BdG) Hamiltonian F^Fo Ui’ 

H = nBdG-^dv 

■HbcIG = {ho + aso(p X a)) -b -b At^, (1) 

where ho = e{p) -b V{r), iFl' = ['i/jj, V'j, V'f) “V't] is the 

Nambu spinor with iptil) being the destruction operator 

for an electron with spin up(down). The kinetic energy 
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term e{p) is given by ^ — /i in a continuum system. 
We consider a 2D wire with p = {Px,Py)- The on-site 
potential is given by V{r), p is the chemical potential 
measured from the bottom of the band, aso is the spin- 
orbit coupling (SOC) strength, B is the Zeeman field and 
A is the proximity-induced s-wave superconducting gap. 
The Pauli matrices cr^ (r^) act on the spin (electron-hole) 
space. 

In the limit of large B, the wire is completely spin po¬ 
larized. Then the low-energy quasiparticles are described 
by an effective p-waye Harnilto nian as discussed in pre¬ 
vious literature.lS^ltSlIllHlIIlIIlll por completeness, we dis¬ 
cuss this limit in Appendi x [P] 

The Hamiltonian in Eq. BJ) is in the Altland-Zirnbauer 
(AZ) symmetry class D {class D for short) in two di¬ 
mension^^ with a topological number Qb, € Z 2 . In 
the absence of SOC along the p-direction, i.e. when the 
®SO Py CTajTz term is set to zero, this Hamiltonian also 
possesses a chiral symmetry, placing it into AZ symme¬ 
try class BDI {class BDI fo r sho rt) with an integer topo¬ 
logical number Qbdi G the thin wire limit, 

i.e. W <C ^sO: chiral symmetry break ing terms are 
0{{W/lso?). Hence, the system in Eq. (IH has an ap¬ 
proximate chiral symmetry. F^F^I^^ ] We show in the next 
section that the class-BDI (chiral) topological number 
Qbdi G Z and the class-D topological number are related 
as Qb, = (—I)‘3bdi (ggg pq^ (^).l^ 


A. Topological index for a disordered multichannel 
s-wave wire 

To obtain the relevant topological index that counts 
the number of the Majorana end states for a RSW wire 
with disorder, we start with the BdG Hamiltonian ^BdG 
in Eq. Q. Following Adagideli et a/.,l^we perform the 
unitary transformation "HbcIG —J' ^BdG ~ ^^^BdG^) 
where U = {\ + i(Jj,){l -b ir^) [1 -b cr^ + (1 - crz)^^] /4. 
Having thus rotated the Hamiltonian to the basis that off- 
diagonalizes its dominant part and leaves the small chiral 
symmetry breaking terms r^crj, in the diagonal block, we 
obtain 


hi'BdG = -Ty ho + aso Px)+Tx {B + A) 

+ T^(JyasoPy- ( 2 ) 

We first set the chiral symmetry breaking term 
Tz o'y OisoPy to zero and focus on A = 0. The eigenvalue 
equation then decouples into the upper and lower spinor 
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components. The solutions are of the form x+ = {4>+i 0)”^ 
and X- = where (j)± obey the following equa¬ 

tion: 


(e(p)cr^ - iPxaso<^x T Bt Acr^) (t>± = 0. (3) 


Here, we have performed an additional rotation (t^ ^ ay, 
ay —>■ —(Tz and premultiplied with ±ax- We note that the 
operator acting on (f)± is not Hermitian. 

We now perform a gauge transformation <p±i^^y) 
(j)±{x,y) with a purely imaginary parameter ina- 
We take Ka to be of first order in aso and identify the 
following terms in the nonhermitian operator in Eq. @ 
in order of increasing power of aso- 


Ho = ho{p; X, y)a^ T ^ T Aua, 


Hi 

H2 


ihtiaPx 

m 



(7z ~\~ hoi^QKct(7x') 


(4) 


where we have indicated the (x, y) dependence of 
ho{p;x,y) through the potential V{x,y). We absorb H 2 
into Hq by redefining p and A. We now identify Ha with 
the inverse of the effective superconducting l ength ^eff , 
setting Ka = = TwasoA/fe with e = y/B'^ — A^. 

With this choice, {Ho,Hi}+ = 0, which allows us to 
write the local solutions as follows: 


{Anfnix,y; e) + Bngn[x,y,e)) 

n 

+ U (-e)e'F'^“^ {Cnfn{x, y, -e) -k Dngn{x, y, -e)), 

(5) 


where (e) are the eigenvectors of the 2x2 matrix ea^ =F 
Actj, with eigenvalue ±|i?| and /„ and gn are the local 
solutions of the equation hgip = e'ljj. The presence of 
a multiple number of local solutions, which is the new 
aspect of the present problem, reflects the multichannel 
nature of the wire. 

We then consider a semi-infinite wire (a; > 0, 0 < 
y < W) described by the Hamiltonian in Eq. Q with 
Gaussian disorder. After going through the steps de¬ 
scribed above, we choose without loss of generality /„ 
to be the decaying and gn the increasing function of x. 
We invoke a well known result from disordered multi¬ 
channel normal state wires and express the asymptotic 
solutions as /„ = e~^’'^Un{x,y) and gn = y) 

where Un{x, y),Vn(x, y) are 0(1) fun ctions as x — » 00 and 
A„ > 0 are the Lyapunov CTponenfs.^^^^^ZIlIIlZIini 

We now focus on a tight-binding system, where the 
number of Lyapunov exponents A^max is finite. (In the 
continuum limit, we have A^max —> 00 .) For the boundary 
conditions at a: = 0, we first extend the hardwall back to 
X = —L' with L' a small value, and consider a normal 
metal in the strip —L' < a; < 0 and Q < y < W (see 
Figure [^(b); in Eq. Q, aso = 0, A = 0, V{x,y) = 0). 
The hardwall boundary condition at a; = —L' can be 
expressed as R- b_^_ = b_ with 6_,_ = (..., A„, Cn, ■ ■ -7, 


b_ = {... ,Bn, Dn, ■■ -7 and R as the extended reflection 
matrix.l^We therefore have 2Nniax boundary conditions, 
leaving 2Nniax of the 4A^i„ax parameters undetermined. 

The boundary conditions at x —)■ 00 require that 4'± 
have only exponentially decaying solutions. We focus on 
the B > A case, yielding real Ka and e. (As discussed in 
ReferenceJi^ ancPSI^ the B < A case yields no solutions.) 
We take Ka > 0 for definiteness. (The following argu¬ 
ments can be extended trivially to the Kq < 0 case.) The 
exponential asymptotic factors in the solutions contain 
a factor of in various sign combinations, affect¬ 

ing the overall convergence at x —)■ 00 . In particular, 
the solutions (j)+ have exponential factors of 
g(«.-eA„(e))x^ g(_„„+A„(-e))a;^ whereas 

the (/)_ solutions have the same form of exponential fac¬ 
tors with the sign of Kq switched. For |kq| smaller than 
all Lyaponov exponents, all and Dn are set to zero 
as they would represent diverging solutions at x —)■ 00 . 
There are therefore 2A^„iax more conditions , bringing 
the total up to 4A^i„ax, to determine a total of dA^max pa¬ 
rameters, yielding only accidental solutions. However, 
for a given n = n*, if min (A„^ (e), (—e)) < Ka < 

max (An, (e), A„, (—e)), there are three growing solutions 
for one of the cj)± sectors and only one for the other sector. 
(If A„, (e) < A„, (—e), the (j)+ sector has the three growing 
solutions and vice versa.) The sector with three growing 
solutions thus has the number of boundary conditions 
increased by one and the other sector has the number 
of boundary conditions decreased by one. If any sector 
has more than 4AAnax boundary conditions in total, then 
there are no solutions for that sector. Therefore, the BDI 
topological number Qbdi S Z is given by the number of 
free parameters, which is equal to 4A^„iax minus the total 
number of equations arising from the boundary condition 
at X = —L'. We obtain: 

Qbdi = X! ® “ ^"(^)) ® (An(-e) - 

n 

- E Q (Cff - An(-e)) 0 (An(e) - Cff ) • (6) 

n 

We see that each Lyapunov exponent pair A„(±e) con¬ 
tributes a topological charge to the overall topo¬ 

logical charge. Hence Qbdi = Qbdo where 

f 4-1 if A„(-e) > > A„(e) 

Qbdi ~ I A„(—e) < < A„(e) 

I 0 otherwise. 

We thus generalize the resuls of Ref.^^ to a multichan¬ 
nel RSW wire. We note, however, that the total num¬ 
ber of Majorana end states for a multichannel RSW 
wire in class BDI, given by IQbdiI, is not equal to sum 
of the Majorana states per Lyapunov exponent pair, 

i-e. IQbdiI ^ IQbdiI- 

We now consider the full Hamiltonian in Eq. Q with 
the chiral symmetry breaking term included. This Hamil¬ 
tonian in two dimensions is in class D and only approxi¬ 
mately in class BDI. The chiral symmetry breaking term 
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pairwise hybridizes the Majorana states described above, 
moving them away from zero energy. However, because 
of the particle-hole symmetry in the topological super¬ 
conductor, any disturbance or any perturbation that is 
higher order in aso can only move the solutions away 
from zero energy eigenvalue in pairs; i.e. for any solu¬ 
tion moving away from zero eigenvalue towards a positive 
value, a matching solution must move to a negative eigen¬ 
value. Therefore, the number of zero eigenvalue solutions 
changes in pairs. Hence, the parity doesn’t change. The 
parity changes, however, every time one of the Lyapunov 
exponents passes through the value of . We therefore 
arrive at the class D topological index Qd = (— 1 )' 3 bdi 

aPl 

Qd = I|sgn(A„(±e)^eff- 1), (7) 

n,± 

indicating that there’s a class D Majorana solution at 
zero energy (Qd = ~1) if there are an odd number of BDI 
Majorana states per edge. Therefore, for the topological 
state of the RSW wire to change from trivial to nontrivial 
or vice versa, it is necessary and sufficient to have Qbdi 
described in Eq. ([^ change by one. The above equation 
thus constitutes the multichannel generalization of Eq.(7) 
of Ref.l^. 

To calculate the topological index Qd in Eq. 0 , we 
relate the Lyapunov exponents in Eq. Q to transport 
properties, namely the mean free path, of a disordered 
wire. We first note that as L —>■ oo, the Lyapunov expo¬ 
nents A„ are self-averaging, with a mean value A„ given 
by 

— Ti 


where is given by [(lT/7ra) arccos (1 — e/2t)\ for 0 < 
e < it and [(W/tto) arccos (1 — (4 — e/2t))\ for it < e < 
8t. The details of the calculation and the dimensionless 
constant are again found in Appendix [a| 

The topological phase boundaries, shown in Eiguresj^ 
and as the bold black lineSjare calculated by equating 
to A„ obtained from Eq. Q and ( [T^ . We thus obtain 
the critical field B* at which the system goes through a 
topological phase transition via thie following implicit 
equation: 


B* = A + 1 (11) 


where /3 = (WaHy-flso?, = ^i± ^{B*)^ + A-^ 

and 




n J 


Equation 0 constitutes the central finding of our pa¬ 
per. It is an analytical expression that determines all 
topological phase boundaries of a multichannel disor¬ 
dered wire. 

An experimentally interesting point is the largest val¬ 
ues of various system parameters that allow a topological 
transition. Using Equations and Q, we estimate the 
upper critical field B*\y , i.e. the minimum value of B 
above which the system is always in a topologically triv¬ 
ial state at a given disorder strength 7 , as 


7max 

5*1^ ~ A (12) 

‘so 


where HeS = ± e, Aj/reg) = [WkF{ties)/T^\, kp = 

n G 1... A(/Xeff) and ^mfp is the MFP of 
the disordered wire.^ We use Eermi’s Golden Rule to 
approximate the mean free path Imfp by calculating the 
lifetime of a momentum state and multiplying it with 
the Fermi speed. We obtain, for a quadratic dispersion 
relation e(p) = I2m — fi, 


I 


-1 _ 

MFP “ 


4to^7 I 
h'^nkp ^ ’ 


(9) 


where is a dimensionless number whose detailed form 
is given in Eq. ( |A5| ). The details of this calculation can 
be found in Appendix [A} 

In order to compare our numerical tight-binding re¬ 
sults with the analytical results obtained through Eq. Q 
and ([^, we also calculate the mean free path 
for a tight-binding (TB) dispersion relation £{kx,n) = 
2t (2 — cos {k^^no) — cos {mra/W)), where t is the hop¬ 
ping parameter, a is the lattice parameter for the TB 
lattice, W is the width of the lattice and kx,n is defined 
through kx^n + ^y,n = with ky^n = wkIW. We obtain 


(&)-' = 


7 


TB\-1 




(Cn) 


( 10 ) 


where = max({A“^}) is the maximum localization 
length achievable in the system. For a hxed nonzero 
disorder, B*\^^q is infinite for a continuum system as 
the localization length increases indefinitely with increas¬ 
ing Fermi energy. For a TB system, the upper critical 
field B*\^^q is finite because the localization length is 
bounded in TB systems. For a clean wire, B*\^^[j is in¬ 
finite for both the TB and the continuum models. 


B. Numerical simulations 

In this section, we obtain the topological index of a 
disordered multichannel wire numerically and compare it 
with our analytical results from the previous section. For 
our numerical simulations, we take the TB form of the 
Hamiltonian in Eq. Q whose details can be found in the 
Appendix [ b 1 We consider a wire of length L ^ ^mfp, or 
Iso, with metallic leads (oso = 0, A = 0 and V{x, y) = 0 
in the leads). We use the results of Fulga et al. to obtain 
the topological quantum number of the disordered multi¬ 
channel wire from the scattering matrices of the wires.^ 
For a semi-infinite wire in the symmetry class D, the 
topological charge is given by Qb = det(r) where r is 
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FIG. 2. (Color online) fi vs. B vs. Qd for a five-channel sys¬ 
tem (compare with Figs. and [^) The background red-white 
colors are obtained using a numerical tight-binding simula¬ 
tion with L = 30000a and W = 5a, while the black lines, 
which represent the topological phase bounda ries, are ob¬ 
tained analytically using Eq. Q. Here, Vo = = 0.2t, 

crso = 0.02/i/ma {Iso ~ 4.08/rm) and A = 0.164t, where 
t = /2ma? and a — O.OlZso is the tight-binding lattice 

spacing. The fragmented nature of the topological phase di¬ 
agram seen in (b) cannot be explained in a p-wave picture. 
See Appendix]^ for a discussion of corresponding experimen¬ 
tal parameters. 



the reflection matrix. For a quasiparticle insulator, this 
determinant can only take the values ±1. However, for 
a finite system this determinant can in general have any 
value in the [—1,1] interval. We obtain the reflection ma¬ 
trix of the TB system in our numerical TB simulations 
using the Kwant librarji^ and then use this relation to 
calculate Qd- We plot the topological phase diagram in 
Figures]^ andwhere the red and white colors represent 
Qd = — 1 and Qd = +1 respectively. 

Figure exemplifies our central result given in 
Eq. pr] ). We find that for a nearly depleted wire 
(Fig. 1^), the topological phase merely shifts to the 
higher values of the chemical potential in agreement with 
Ref.l^. For higher chemical potentials/doping, we ob¬ 
serve a fragmented topological phase diagram (Fig. Ep)- 
We find good agreement with our analytical results from 
Eq. (111. We note in passing that, this fragmenta¬ 
tion cannot be explained by a simple p-wave picture 
as these topological phases arise despite the incomplete 
spin-polarization of the wire under a low magnetic field. 
For a full phase diagram over the entire bandwidth, but 
for slightly different material parameters, see Figure [^ 
where the reentrant phases are apparent. 

In Fig.[^ we plot the topological n umber Qd as a func¬ 
tion of /r and the disorder strength for a constant 

Bzeeman Over the full TB bandwidth. The reentrant na¬ 
ture of the topological phase diagram can also be seen in 
this plot, for example, by following the /r = 1.5 line as 7 
is increased. As the disorder strength increases, series of 
topological transitions occur, similar to the PW wire.l^ 
However, unlike the PW wire, the number of transitions 
is given hy N{fi + e) + N{^ — e) rather than fV(/r), with 
Nip) defined as N{iies) = [lTfc_F(/reff)/ 7 rJ. For further 
discussion of the emergence of effective p-wave picture at 
high magnetic fields, see Appendix [C| 
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FIG. 3. (Color online) p vs. Vo = •Jyja? vs. Q for a 
multichannel RSW wire. The black lines, which represent 
topological phase boundaries, are obtained analytically us¬ 
ing Eq. 0. The background red-white colors are obtained 
using tight-binding numerical simulations with L = 60000a. 
In both cases, W = 4a, aso = 0.0151i/ma, A = 0.20t and 
B = 0.35t, where t = T? j^maP is the tight-binding hopping 
parameter and a is the TB lattice spacing. See Appendix [B| 
for a discussion of corresponding experimental parameters. 


III. CONCLUSION 

In summary, we investigate the effect of disorder in 
multichannel Rashba SOC proximity-induced topological 
superconductor nanowires (RSW nanowires) at experi¬ 
mentally relevant parameter ranges. We derive formulae 
that determine all topological phase boundaries of a mul¬ 
tichannel disordered RSW wire. We test these formulae 
with numerical tight-binding simulations at experimen¬ 
tally relevant parameter ranges and And good agreement 
without any fitting parameters. We show that there are 
additional topological transitions for the RSW wires lead¬ 
ing to a richer phase diagram with further fragmentaliza- 
tion beyond that of the p-wave models. 
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Appendix A: Mean free path 

We consider a long wire along the x-axis, having a 
length of L along the x-direction and a width of W along 
the 2 /-direction and metallic leads at the end, with a 
Gaussian disorder of the form (^(r) C(r')) = 7(5(r — r'). 
We obtain the ensemble average of the matrix element be¬ 
tween the and /*** transverse channels as k{kx,n) —)■ 
k' as 



FIG. 4. + 1) vs. N. 



We then use Fermi’s Golden Rule to calculate the inverse 
lifetime of a momentum state k, : 

/,-i \ _ ^ 27r 7 

g,.J X 

(l + %^) P(£.')- (A2) 


value of (n just below the transition N —>• A^-|-l (denoted 
(n^n+i) is plotted in Figure]^ 

We now derive the MFP for a TB dispersion relation 
given by 

£{kx,n) = (2 — cos {kx^no) — cos {mra/W)). (A6) 

The number of channels is given by iV = 
[(IF/7ra) arccos (1 — e/2f)J for 0 < e < 4t and by iV = 
[(IF/7ra) arccos (1 — (4 — e/2t))J for 4t < e < 8t. The 
resulting disorder-averaged inverse MFP reads: 


where Sk gives the dispersion relation and p{sk) is the 
density of states. We then sum over the initial and final 
states k' in Eq. ( |A2[ ) to obtain the total inverse MFP: 

(^MFp) = ^ (^MFP(fc^,n^fei.o) (A3) 


We first apply Eq. (A3) to a free electron dispersion of 
the form s{k) = h^Wj2m = h‘^/2m{k'^^^ + n^7r^/W^) 
for n G l,...,iV where N{p,ee) = [WkF{e)/TT\. The 
resulting total ensemble-averaged inverse MFP is 


/,_! \ _ 27IF ^^ ^ 6 nl ^ 

{^MFP) -1^2^ I + T j 

n=l ;=i •> ' V / 


TT 



± sj2meir9 - jW^) 
\j2me/k? — rP-TT^ jW^ \/2m£/h? — jW"^ 


Arn?"f 

ti^TTkF 


CAr\ 


(A4) 


((^MFP) (C™) ^ (A7) 

where the dimensionless is given by 

n —1 
N N 

(AS) 

n —1 />n 


Here, r]™ = | sin a)| ^ 
using Eq. (A6) 


and sin {kx,n) is obtained 


Appendix B: Numerical tight-binding simulations 

We start W obtaining the TB form of the RSW BdG 
HamiltoniarpS in Eq. Q in the usual way using finite 
differences (see, for example. Ref. It reads: 


where kp = yj2/m£2fW is the Fermi wavevector, 

^ N _ N N 

cy = ^E^" + 2iV^^77„77i, (A5) 

n —1 n —1 l^n 

_ 1 

and r]n = ~ > ic^ agreement with Eq.(8) in 

the supporting online material of Rieder et al^^. The 


■^EdG = [(4^ + V(x, y) - y{x, y)) + Hz 


^{x,y)Tx\ \x,y) {x,y\ 
i . 

- -^aso[x,y)Tzay 


+ -aso{x,y)Tzax 


+ h.c. 


\x + a,y) {x,y\ 

\x,y + a) {x,y\ 
(Bl) 
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where t = /2ma?' is the hopping parameter, V(x,y) 

is the Gaussian random potential, y) is the rele¬ 
vant gate potential, i?z is the Zeeman field, A{x,y) is 
the s-wave superconducting pairing (taken to be real), 
<^So(x,y) is the effective Rashba SOC due to proxim¬ 
ity effect and a is the lattice constant for the TB lattice. 
Here, V{x,y), Bz, A(x,y) and asoix^y) are nonzero only 
within the scattering region. Bz, asoix,y) and A(x,y) 
are constant within the scattering region except for the 
values of aso(x,y) in the scattering region-lead bound¬ 
ary, where we take it to be half of its value in the bulk. 

The experimental values for InSb nanowires quoted 
in Mourik et al^ are aso = 0.2 eVA, Iso 2000A, 
A = 0.25meV, Ez/B = 1.5meV/T, to* = 0.015TOe and 
agQm^l2h‘^ ~ O.OdTOeV. We employ these values verba¬ 
tim, except for Iso (and correspondingly, aso), for which 
we use parameters much more accessible experimentally. 

We use the Kwant librarjJ^ to obtain the topologi¬ 
cal phase diagram in our numerical plots. The Kwant 
library can extract the scattering matrix (S'-matrix),!^ 
and therefore the reflection matrix (r-matrix) for a given 
tight-binding system with leads. The topological in¬ 
dex Qo can be obtained from the r-matrix through 
Qo = det(r) (see Ref.^. 

The numerical parameters quoted in the caption of 
Fig.i correspond to f = l.STOeV, a = 40.8nm, Iso = 
4.08^m and a = 6.3 x 10“® c. Disregarding screening, 
a Zeeman Energy of, say, 0.35f on the plot would corre¬ 
spond to a magnetic field 0.35T, a value easily accessible 
by the experiment. In Figures[^[^[^and|^, Iso = 6.0/im, 
t = O.TmeV, a = 60.0nm and a = 4.2 x 10“® c. A Zeeman 
energy of 0.35t corresponds to B = 0.17 T. 

The TB form of the effective PW Hamiltonian of 
Eq. (Dl) used in Appendix]^ is as follows: 


/TB 


[it + V{x,y) - y,{x,y)] 

\x,y) {x,y\ 

-f 

-fT/; - ^Aeg{x,y) Tx 

\x + a,y) {x,y\ 

-f 

-tT^ - ]^Aeg{x,y)Ty 

\x,y + a) {x,y\ 

-1- h.c. 

( 


(B2) 


We use numerical values similar to the RSW case 
in our PW simulations, except to impose Aes = 
Aaso/Vi?2-A2. 


Appendix C: Topological phase diagram over the full 
bandwidth 


photonic metamaterial^^ simulating topological proper- 
ties.l ^ Al l analytical plots are produced using Eq. Q 
(Eq. ( |D2[ ) for the PW case), but using a TB dispersion 
relation for e{p) in the relevant expressions. All of the 
numerical results are obtained using a TB simulation uti¬ 
lizing Kwant software, as discussed in the main text. 


Figure depicts the analytically calculated topolog¬ 
ical phase diagram for an RSW wire as a function of 
/i and the disorder strength, for various magnetic field 
strengths. The transition between a RSW wire and a pair 
of oppositely polarized PW wires can be seen as increas¬ 
ing magnetic field polarizes the system. The topological 
order is less robust against disorder for higher magnetic 
fields, because the coherence length becomes longer with 
increasing B. This is the reason why the spin polarized 
regimes where PW model applies is typically less robust 
than the lower field regimes where both spin species exist 
as seen in Fig. [^a) and[^c) or (d). In order to complete 
the discussion, we also present an analytical plot (Fig¬ 
ure 1^ for an RSW wire for which B is greater than the 
subband spacing but less than the bandwidth. While this 
regime is experimentally very hard to achieve, it is useful 
for comparing the PW and the RSW regimes. The verti¬ 
cal blue line denotes the bottom of the higher energy spin 
band beyond which both spin species exist. We note that 
the critical disorder strength increases with the chemical 
potential, hence spin-polarized regime, which appear at 
lower chemical potential values, is less robust against dis¬ 
order. 


In Figure the analytically calculated phase diagram 
of a wire with W = 4a is plotted with increasing disor¬ 
der. We see that the phase diagram gets fragmented as 
number of channels are increased. We also note that for 
a given amount of disorder, there is a maximum Zeeman 
field Hmax above which no topological order is present. 
The reason is that in our numerical TB simulations, the 
localization length is not a monotonous function of en¬ 
ergy. It grows (with increasing energy) until the middle of 
the band, and after that it decreases as the energy comes 
closer to the band edge. This places an upper magnetic 
field limit to topological regions since the superconduct¬ 
ing coherence length monotonically increases with B. For 
a pure quadratic dispersion, the upper limit is given by 
the limitations of the approximations of Fermi’s Golden 
Rule and would increase indefinitely with increasing en¬ 
ergy as discussed in the main text. We note that the up¬ 
per limit discussed here has a different origin than that 
discussed by Ref.^ for finite-length wires. 


In this section, we present plots of the topological 
phase diagram that we obtain analytically from Eq. 0 
using a TB dispersion relation (see Section 0 over the 
full bandwidth. Although only the low y regions in our 
plots correspond to experimentally relevant nanowires, 
the full bandwidth range would be important for sys¬ 
tems that are inherently TB, such as atomic chains^ or 


We finally present the full TB bandwidth version of 
Fig.H with slightly different material properties, here 
in Fig. 1^ This figure is the numerical simulation result 
that matches the last of the analytical plots in Fig. 
The relevant numerical values are given in each of the 
Figures’ captions. 








1.2 






Gate potential (^/t) 


I—I ^ 


0 


FIG. 5. (Color online) /r vs. Vo ~ y/'ila? vs. Qd for a 
multichanneled RSW wire for different B, obtained analyti¬ 
cally using Eq. Q. a), b) Low magnetic field B > A limit 
requires a full RSW model and topological order can survive 
up to high disorder strengths, c), d) The spin-polarized sys¬ 
tem can be described by a PW model and topological order 
is completely destroyed with less disorder. Here, W = 4a, 
aso = 0.015fi/ma and A = 0.20t where t = T? 1‘lma} and 
a is the tight-binding lattice spacing. See Appendix [B| for a 
discussion of corresponding experimental parameters. 


Appendix D: Topological phase diagram for 
multichannel effective p-wave nanowires with 
disorder 


In this Appendix section, we present the effects of dis¬ 
order on PW wire s, which is a system previously stud¬ 
ied in literature, I 22 H 2 IIIIHIIIIZIIII for completeness and for 
comparison with the results of our paper for disordered 
multichannel RSW nanowires. We start with the Hamil¬ 


tonian in Eq. (Dl) and present the topological charge 


in Eq. (D2). We plot the topological phase diagram for 


a PW wire as a function of p and disorder strength for 
a fixed Ll 2 eeman E and compare this plot with its 
analogue for RSW wires (Fig. |^. 

The BdG Hamiltonian for an effective p-wave wire with 
spatially homogeneous effective SOC strength is 


u 

Jl 


PW 

BdG 


= e{p)Tz -b AeffP • r. 


(Dl) 


Note that Aes has units of velocity while A in Eq. Q 
has units of energy. This effective SOC strength is re¬ 
lated to the corr esponding RSW superconducting gap by 
Aoff = Aaso /— A^.li^ We consider a Gaussian dis¬ 
order of the form (H(r)H(r')) = j6{r — r') for r,r' in 
the wire, with 7 as the disorder strength and (H(r)) = 0 . 
This Hamiltonian is useful for comparison with the fully 
polarized limit of the RSW case . 

The Hamiltonian in Eq. (Dl) is in Altland-Zirnbauer 
(AZ) symmetry class D in two dimension^^ with a Z 2 
topological number. This Hamiltonian also possesses a 



FIG. 6. (Color online) p vs. Vo = y/'ila? vs. Qd for a 
multichanneled wide RSW wire, obtained obtained analyti¬ 
cally using Eq. 0- with W = 77a. Here, aso = 0.015fi/ma 
(Iso = 100a), A = 0.20t and B — 0.205t with the hopping 
parameter t = f2ma^ = 0.7 meV and the lattice spacing 
a = 60 nm. The blue vertical line at /r = e = \/B'^ — A^ is 
the bottom of the second spin band. 


If f ff f f f 

a) = O.Of 


, 


ff f 



b) = O.U 


“I-1-1-1-i-1-r 




d) = 0.3t 

I I I I I 



c) ■s/'ya? = 0.2< 
I’ll 


1234567 1234567 

Gate poteatial (^4/i) 


FIG. 7. (Color online) /i vs. B vs. Qd for varying disorder 
strengths for an RSW TS with Gaussian disorder, analytically 
calculated using Eq. ([^ for a four-channel TB system. Sub¬ 
figure c) matches the numerical data shown in Fig. . The 
parameters used are aso = 0.0151i/ma and A = 0.2f, where 
t = 12ma^ a is the lattice spacing. See Appendix |b| for a 
discussion of corresponding experimental parameters. 


chiral symmetry, broken by the A^epyTy term. If this 
ter m is set t o zero, the Hamiltonian is also in class 
Bd 14 ^E 6 | 58 | 59 | Slaving a Z topological number. (ID wires 
trivially satisfy this condition.) In the thin wire limit, 
i.e. Aeff <?; H/mW, the chiral symmetry breaking term 
is O ((mAefflT/fi.)^). The wire in class BDI can have an 
integer number of Majorana fermions at its ends. The 
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0 I-1-1-1-1-1-1-!-i H _1 

012345678 
Gate potential {n/t) 


chiral symmetry breaking term pairwise hybridizes these 
solutions. Hence the chiral topological number Qbdi G Z 
and the class-D topological number Qd G ^2 are related 
as Qd = 


In order to solve the Schrodinger equation H 
at i? = 0 to obtain the Lyapunov exponents, we follow 
Adagideli et al¥^ to off-diagonalize the Hamiltonian and 
apply an imaginary gauge transformation. This allows 
us to re-express Qbdi in terms of A„!^ 


FIG. 8. (Color online) p vs. B vs. Qd for a four-channel sys¬ 
tem (compare with Figs. I^andj^. The black lines, which rep¬ 
resent topological boundaries, are obtained analytically using 
Eq. 0. The background red-white colors are obtained us¬ 
ing tight-binding numerical simulations. The parameters are 
Vo = 0.2t, A = 0.2t and aso = 0.015h/ma. See Appendix [b| 
for a discussion of corresponding experimental parameters. 


N 

Qbdi = ^ 0 

n—1 



(D2) 



FIG. 9. (Color online) /r vs. yj'y/a? vs. Q for a multichan- 
neled PW wire with dimensions W = Aa and L — 60000a 
{L used only in the numerical tight-binding code) and with 
Qso = O.Olh/ma, where a is the tight-binding lattice spacing. 
The red-white colors in the background are obtained numer¬ 
ically with a tight-binding meth od wherea s th e black solid 
lines are obtained using Eq. (D2l with Eq. (lOl. 


where N = [W/tt^J 2rr^^/W\ and \x\ is the usual floor 
function. We obtain A„ again using Eq. (0. We obtain 
^MFP using Fermi’s Golden Rule (see Appendix [a]) first 
for a quadratic dispersion relation and then for a TB 
dispersion relation. 


We compare the results found using Eq. |D2| with those 
obtained by numerical simulations in Figure!^ and find an 
excellent fit over the whole TB bandwidth. In a clean PW 
wire (\A)7^ = 0), Majorana modes appear if N is odd 
and Majorana states fuse to form ordinary Dirac fermions 
if N is even. This behavior survives up to a finite disor¬ 
der strength (see Fig. 0. As in the case of RSW wires, 
further increase of the disorder strength gives a series 
of transitions between non-trivial and trivial topological 
phases as each A„ increases and crosses While both 
multichanneled RSW and PW wires feature reentrant be¬ 
havior, we see that there are additional transitions for the 
RSW wires leading to a richer phase diagram (compare 
Figures]^ an di, in agreement with our analytical results 
presented in Eq. 0- 
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